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The half-life for double beta decay of 150 Nd has been measured by the NEMO-3 experiment at 
the Modane Underground Laboratory. Using 924.7 days of data recorded with 36.55 g of 150 Nd the 
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half-life for 2^/3/3 decay is measured to be T 2 J 2 = (9.1llg;g(stat.) ± 0.63 (syst.)) x 10 18 years. The 
observed limit on the half- life for neutrinoless double beta decay is found to be T°J 2 > 1-8 x 10 22 years 
at 90% Confidence Level. This translates into a limit on the effective Majorana neutrino mass of 
(m v ) < 4.0 — 6.3 eV if the nuclear deformation is taken into account. We also set limits on models 
involving Majoron emission, right-handed currents and transitions to excited states. 

PACS numbers: 14.60.Pq, 23.40.-s 



Experiments studying atmospheric, solar, reactor and 
accelerator neutrinos have established the existence of 
neutrino oscillations as a direct evidence for a non-zero 
neutrino mass. These experiments, however, cannot dis- 
tinguish between Dirac or Majorana neutrinos. They 
also provide no information on the absolute neutrino 
mass scale, since oscillations experiments measure the 
square of the mass difference between neutrino states. 
The half-life of neutrinoless double beta decay (OV/3/3) is 
inversely proportional to the effective Majorana neutrino 
mass squared, (m u ) 2 . Observation of this process would 



therefore directly constrain the neutrino mass scale and 
would be unambiguous evidence for the Majorana nature 
of neutrinos. The 0V/3/3 process also violates lepton num- 
ber and is therefore a direct probe for physics beyond the 
standard model of particle physics. 

The search for neutrinoless double beta decay of 
neodymium-150 ( 150 Nd) using the NEMO-3 detector is 
of special interest since 150 Nd is a candidate isotope 
for SuperNEMO a next generation double beta de- 
cay experiment based on the NEMO-3 concept, and the 
SNO++ experiment at SNOLAB 0. Its main advan- 
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tages are the high Qpp value for double beta decay, 
Q/3/3 = 3.368 MeV, which lies above the typical ener- 
gies for many background sources, and the large phase 
space factor. The 2v(3(3 half-life of 150 Nd has previously 
been measured using a Time Projection Chamber 

The NEMO-3 experiment has been taking data since 
2003 in the Modane Underground Laboratory (LSM) lo- 
cated in the Frejus tunnel at a depth of 4800 m water 
equivalent. The detector has a cylindrical shape with 
20 sectors that contain different isotopes in the form of 
thin foils with a total surface of about 20 m 2 Q. In 
addition to -7 kg of 100 Mo and -1 kg of 82 Se, the de- 
tector contains smaller amounts of other isotopes. The 
neodymium source foil is composed of Nd2 03 with an 
enrichement of (91 ± 0.5)%, corresponding to a 150 Nd 
mass of 36.55 ± 0.10 g. On each side of the foils is a 
^50 cm wide tracking volume comprising a total of 6180 
drift cells operated in Geiger mode with helium as drift 
gas. A 25 Gauss magnetic field created by a solenoid 
provides charge identification. The calorimeter consists 
of 1940 plastic scintillators coupled to low radioactivity 
photomultipliers. For 1 MeV electrons the energy res- 
olution (FWHM) ranges from 14.1% to 17.7% and the 
timing resolution is 250 ps. A cylindrical coordinate sys- 
tem (r, z, <fi) is used with the z axis pointing upwards. 

The data used in this Letter have been recorded be- 
tween February 2003 and December 2006, correspond- 
ing to 924.7 days of data taking. Signal and background 
events are generated with the GENBB generator and simu- 
lated using a GEANT-based simulation @ of the detector. 
All Monte Carlo events are processed by the same recon- 
struction programs as the data. A detailed description 
of the analysis can be found in Q • 

The 2v(3(3 events are expected to have two electrons in 
the final state and they are therefore selected by requir- 
ing two tracks with a curvature consistent with a negative 
charge. Each track has to be matched to a separate en- 
ergy deposit in the calorimeter greater than 0.2 MeV. The 
z component of the distance between the intersections of 
each track with the plane of the foil should be less than 
4 cm and the transverse component less than 2 cm. Both 
tracks must originate from the first layer of the Geiger 
cells. The scintillator measurement of the time-of-flight 
(TOF) for both electron candidates must be consistent 
with the hypothesis that the event originates from the 
source foil. After this selection 2789 events remain. 

The background sources are divided into two cate- 
gories, depending on their origin . Background events 
originating from the radioactive impurities in the source 
foils are called internal background. (3 emitters can pro- 
duce /3/3-like events through three mechanisms: (i) (3 de- 
cay accompanied by electron conversion; (ii) M0ller scat- 
tering of a (3 decay electron, and (iii) (3 decay to an ex- 
cited state followed by Compton scattering of the de- 
excitation photon. Another internal background source 
is 207 Bi, most likely present due to a contamination of the 



150 Nd source foil during production. This isotope decays 
to excited 207 Pb through electron capture and 207 Pb sub- 
sequently emits two electrons via electron conversion. 

The second source of background is called external and 
is caused by electrons or photons anywhere outside the 
source foils. The main source of external background is 
radon. Radon decays to 214 Bi via two a decays and one 
(3~ decay. The decay of 214 Bi to 214 Po is generally ac- 
companied by one electron and several photons, which 
can mimic a (3(3 event through conversions. A radon pu- 
rification facility was installed about half-way through 
the data taking period presented in this Letter. This re- 
duces the radon induced background for this analysis by 
about a factor of six. Other external background contri- 
butions are found to be small. 

The background activities are determined by measur- 
ing control decay channels that are independent of the 
signal topology using a full MC simulation of the back- 
ground processes and the detector. A summary of the 
measured background activities are given in Table [U 

In the ej channel, 207 Pb from the 207 Bi decays via the 
emission of an electron and a photon from the strongly- 
converted energy transition in excited 207 Pb. 152 Eu de- 
cays into excited 152 Gd through (3 decay which de-excites 
into the ground state via photon emission. The e"f events 
are selected by requiring that exactly one negatively 
charged particle with a track length greater than 50 cm is 
found in the sector containing the 150 Nd foil. The track 
must be associated with an isolated scintillator hit with 
energy greater than 0.2 MeV. It must originate from the 
150 Nd foil and must have a hit in the first layer of the 
Geiger cells. The scintillator measurement of the time- 
of-flight of the electron candidates has to be consistent 
with the hypothesis that the decay occurred in the source 
foil. The photon is identified by requiring that there is a 
second scintillator hit with no associated track and with 
energy greater than 0.2 MeV. The energy sum of the all 
other clusters not associated to the electron or the photon 
must be less than 0.15 MeV. The opening angle between 
the electron and the photon is required to be cos 9 < 0.9. 

To obtain the background normalisation for 207 Bi and 
152 Eu, the normalization of these background contribu- 
tions is fitted to spectra of the electron energies, E e , and 
the photon energies, E 1 , shown in Figure [TJ The contri- 
butions from other background sources are fixed in the fit. 
The 228 Ac contribution is taken from the 208 T1 compo- 
nent in the high energy tail of the E 1 distribution, since 
208 T1 originates from 228 Ac decays. The 212 Bi compo- 
nent is normalised with respect to the measurement of its 
decay product 208 T1. The external, the radon induced, 
214 Bi and 214 Pb backgrounds are set to their indepen- 
dently measured values. 

The isotopes 234m Pa and 40 K undergo (3 decay. Their 
activities are therefore measured with single electron 
events. The same selection for the electron is applied 
as in the e-y channel, apart from an additional require- 
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FIG. 1: (color online) Energy distributions of the a) electrons 
and b) photons in cy events. The data with statistical uncer- 
tainties are compared to the sum of the expected background. 



ment E e > 0.5 MeV. The 234m Pa and 40 K contributions 
are obtained from a fit of the normalization to the E e 
distributions. The external background is set to the val- 
ues obtained from independent measurements. The 207 Bi 
component has been measured in the ej channel. The E e 
distribution with the fitted 234rra Pa and 40 K contributions 
is shown in Figure The distribution of the event ver- 
tices for the single electron events, calculated as the elec- 
tron track's intersection with the source foil, shows small 
regions of high activity consistent with 234m Pa contami- 
nation which are removed from the analysis. The activ- 
ities given in Table [U are measured after the removal of 
these regions. 

Due to the resolution of the tracking detector, events 
from neighboring source foils ( 48 Ca, 96 Zr, 100 Mo) can 
be reconstructed as originating from the 150 Nd foil. The 
total number of background events from neighboring foils 
is estimated to be 126 ± 12 from MC simulations, where 
the uncertainty is dominated by the uncertainty on the 
half-life of these isotopes. 

The distributions of the energy sum of the two elec- 



FIG. 2: (color online) Distribution of the electron energy in 
single electron events. The data are compared to the sum of 
the expected background. 



Background 


A (mBq/kg) 


N bB 


N bg for 
£ sum > 2.5 MeV 


ib2 Eu 


54 ±6 


19 ±1 





154 Eu 


22 ±2 


9± 1 





208 rpj 


10 ±2 


46 ±4 


3.5 ±0.9 


228 Ac 


27.8 ±5.6 


52 ±5 





212 Bi 


27.8 ±5.6 


32 ±4 





207 Bi 


231 ± 10 


138 ±6 





214 Bi 


3.3 ±0.8 


13 ±3 


0.6 ±0.2 


214p b 


3.3 ±0.8 


6± 1 





40 K 


213 ±10 


66 ±6 





234 m p a 


47 ±2 


143 ±6 





External and 








radon induced 




53 ±11 


4.8 ±0.8 


48 Ca, 96 Zr, 100 Mo 




168 ± 19 


0.12 ±0.02 


sum 




746 ± 30 


9.0 ± 1.2 


data 




2789 


29 



TABLE I: Summary of the measured background activities, 
the expected number of background events and the observed 
number of events for the selected data set and in the high 
energy region, -B aU m > 2.5 MeV. 



trons and the opening angle between them are shown in 
Figure [3] The data are in good agreement with the sum 
of the background and the 2vj3(3 signal distributions. We 
therefore use the data sample to measure the 2v(3j3 half- 
life, TfY 2 , of 150 Nd. The efficiency of the 2vf3/3 event 
selection is 7.2%. After background subtraction we ob- 
tain 
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(9.11^;^(stat.) ± 0.63 (syst.)) x 10 iS y . (1) 



The systematic uncertainty on the sum of the internal 
and external background is 4.3% which translates into an 
uncertainty on the 2z//3/3 half-life of 1.6%. This includes 
the uncertainty from the background fits and from the 
measurement of the activity of 208 Tl using two different 
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FIG. 3: (color online) Distributions of a) the energy sum of 
the two electrons, E suln , and b) the angle between the two 
electrons, cos#, for data compared to the sum of the back- 
ground and 2vf3f3 signal expectations. 



populated using Poisson simulations of the two hypothe- 
ses. Systematic uncertainties are treated as uncertain- 
ties on the expected numbers of events and are folded 
into the signal and background expectations via a Gaus- 
sian distribution. Correlations between systematic un- 
certainties are taken into account. The value of the con- 
fidence level, CL S , is denned as CL S = CL s+ b/CLb, 
where CL s+ t, and CLb are the confidence levels in the 
signal-plus-background and background only hypotheses, 
respectively. The limits are calculated by scaling the sig- 
nal until 1 — CL S reaches 0.9. 
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FIG. 4: (color online) Distribution of the energy sum of the 
two electrons, -E sum , for E suln > 2.5 MeV. The data are com- 
pared to the total background, consisting of internal and ex- 
ternal background and the 2i//3/3 expectation. A MC simu- 
lation of a Ov/3/3 signal with a half-life of 1.8 x 10 21 years, 
corresponding to ten times the number of events expected for 
the observed 90% C.L. limit, is also shown. 



decay channels. The uncertainty on the tracking effi- 
ciency is 5.7%. Varying the TOF requirement leads to 
a 1% uncertainty. An uncertainty of 3% is due to the 
uncertainty on the position of the 150 Nd foil in the de- 
tector. The uncertainty on the energy measurement was 
studied by smearing the individual electron energies by 
2% which yields a systematic uncertainty of 1.5% on the 
2v(3p half-life. 

Since no significant excess is observed in the E aVLra dis- 
tribution, a limit is set on the half-life for neutrinoless 
double beta decay, T^Y 2 , using the CL S method Only 
the shapes and not the normalization of the full E sum 
distribution are used to discriminate signal and 2v(3(3 
background. The normalization of the -E sum distribu- 
tion for the other backgrounds is fixed to the value given 
in Table |U All limits in this Letter are calculated by 
utilizing a likelihood-fitter [1(3] that uses a log-likelihood 
ratio (LLR) test statistic method. Two hypotheses are 
defined, the signal-plus-background hypothesis and the 
background-only hypothesis. The LLR distributions are 



The total efficiency for 0v/3f3 events after applying all 
selections is (19±1)%. The uncertainties on the efficiency 
of the signal and the background are assumed to be fully 
correlated. The £ sum distribution for E sum > 2.5 MeV 
is shown in Figure [4] for data compared to the total 
background, which consists of internal and external back- 
ground and the 2i/f3/3 expectation. A MC simulation of 
a OV/3/3 signal is also shown. Signal and background are 
well separated, demonstrating the advantages of using 
150 Nd for future 0^/3/3 searches. The observed limit on 
the half-life is T°J 2 > 1 .8 x 10 22 years at 90% Confidence 
Level (CL). It is consistent with the median expected 
limit. 

This limit on the half-life is converted into a limit on 
the effective Majorana neutrino mass, (m„), using nu- 
clear matrix elements (NME). In the currently available 
QRPA-like calculations spherical symmetry of the nu- 
cleus has been assumed. If the upper and lower limits 
on the calculated NME, which also include the uncer- 
tainties in the weak coupling constant <ju, are taken into 
account flj , the experimental lower limit on the half-life 
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of 150 Nd translates into a limit (m v ) < 1.5 - 2.5 eV. 

Taking into account the nuclear deformation will mod- 
ify this conclusion. The suppression of the NME for 
150 Nd has been estimated to be a factor 2.7 in the 
case of nuclear deformations derived from laboratory mo 
This increases the upper limit to (m 
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< 



ments 

4.0 — 6.8 eV which is consistent with the limit derived 
using the NME of a pseudo SU(3) model and the 
PHFB model [3], which also include the effect of nu- 
clear deformation. Further progress in the calculation of 
the NME for (V/3/3 decay of 150 Nd is therefore urgently 
required. 

Neutrinoless double beta decay can also proceed 
through the emission of one or two Majorons (x°) [Hj]. 
These decays are characterized by the spectral index n 
which leads to a modification of the energy sum distri- 
bution by a factor (Qpp — E sum Y\ We consider models 
with n — 1, 3, 7 (l5| and n — 2 [16j |. Limits on n = 3, 7 
models are less sensitive than n = 1, 2 since the shape of 
the -E sum distribution is similar to the one for 2v[3(3 decay 
(n = 5). The CL S method is applied to the E sum dis- 
tribution in the same way as described previously. Since 
the normalisation of the 2v(3f3 background can not be in- 
dependently determined, it is left unconstrained in the 
fit determining the limits. 

We also set limits for Qv[3f3 processes involving right- 
handed currents (V+A). In addition, the 0vf3(3 decay 
can also proceed through various excited states (0^ or 
2^). Limits on all these modes of neutrinoless double 
beta decay are shown in Tabic [TTJ They cither improve 
previously published limits [4j, ll7j significantly or repre- 
sent the first limits for this isotope (n = 2,3,7). The 
limit on n — 1 Majoron emission can be be translated 
into a limit on the neutrino-Majoron coupling of g ee < 
(0.64 - 1.05) x 10~ 4 without and g ee < (1.7 - 3.0) x 10" 4 
with nuclear deformation. These limits are comparable 
with the limits obtained for 82 Se and 100 Mo with about 9 
and 70 times the exposure, defined by isotope mass times 
observation time 181. 



surement of the half-life of double beta decay of 5U Nd 
to date, yielding a value of Tfj 2 = (9.11^^22 (stat.) ± 
0.63 (syst.)) x 10 18 years. This value is slightly more 
than two standard deviations higher than the previously 





Qvf3f3 


Majorons 


0+ «m„)) 


0+(V+A) 


2f 


0r + 


n = 1 


n = 2 


n = 3 


n = 7 


limit 


18.0 


10.7 


2.4 


0.24 


1.52 


0.54 


0.22 


0.047 



TABLE II: 90% C.L. limits on the half-life, T 1/2 , in units of 
10 21 y, for different modes of neutrinoless double beta decay. 



In summary, we have presented the most precise mea- 



(6.75+^2(stat.)±0.68 (syst.)) 



measured value T^Y 2 
10 18 years We have significantly improved limits on 
the half- life of different modes of neutrinoless double beta 
decay for 150 Nd and we have provided the first limits for 
several Majoron models using 150 Nd. Our measurements 
demonstrate that 150 Nd is an excellent candidate isotope 
for future double beta decay experiments. 
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